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Abstract. Like-charged macroions attract each other as a result of strong electrostatic correlations in the 
presence of multivalent counterions or at low temperatures. We investigate the effective electrostatic in- 
teraction between i) two like-charged rods and ii) two like-charged spheres using the recently introduced 
strong-coupling theory, which becomes asymptotically exact in the limit of large coupling parameter (i. e. for 
large counterion valency, low temperature, or high surface charge density on macroions) . In contrast to pre- 
vious applications of the strong-coupling theory, we deal with curved surfaces and an additional parameter, 
referred to as Manning parameter, is introduced, which measures the ratio between the radius of curvature 
of macroions to the Gouy-Chapman length. This parameter, together with the size of the confining box en- 
closing the two macroions and their neutralizing counterions, controls the counterion-condensation process 
that directly affects the effective interactions. For sufficiently large Manning parameters (weakly-curved 
surfaces) , we find a strong long-ranged attraction between two macroions that form a closely-packed bound 
state with small surface-to-surface separation of the order of the counterion diameter in agreement with 
recent simulations results. For small Manning parameters (highly-curved surfaces) , on the other hand, the 
equilibrium separation increases and the macroions unbind from each other as the confinement volume 
increases to infinity. This occurs via a continuous universal unbinding transition for two charged rods at 
a threshold Manning parameter of £ c = 2/3, while the transition is strongly discontinuous for spheres be- 
cause of a pronounced potential barrier at intermediate distances. Unlike the cylindrical case, the attractive 
forces between spheres disappear for increasing confinement volume due to the complete de-condensation 
of counterions. Scaling arguments suggest that for moderate values of coupling parameter, strong-coupling 
predictions remain valid for sufficiently small surface-to-surface separations. 
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1 Introduction 

Electrostatic interactions play a prominent role in many 
soft-matter and biological systems, since a huge class of 
macromolecules are water-soluble, and thus bear electric 
charges in aqueous solutions. These macroions may be 
polyelectrolytes such as DNA, or charged colloidal parti- 
cles, that are surrounded by their neutralizing (oppositely- 
charged) counterions. In recent years, mounting evidence 
from both experiments [fj2 ,3,4,5, 6j and numerical simu- 
lations [7TiaiiMTT)in^^ 

I24U25| showed that like-charged macroions can attract each 
other via effective forces of electrostatic origin, notori- 
ous examples of which are formation of dense packages 
of DNA molecules (DNA condensates) and large aggre- 
gates of charged colloidal particles 18 , 19 , 20 . These obser- 
vations indicate the appearance of like-charge attraction 



in strongly-charged systems, i.e. when multivalent coun- 
terions are present, macroions are highly charged, or the 
strength of electrostatic interactions is enhanced by main- 
taining the system at low temperatures or in a medium of 
low dielectric constant 



On the other hand, there has been a great deal of 
theoretical studies aimed at revealing possible underlying 
mechanisms, which can lead to such counter-intuitive in- 
teractions i7iisii9TrLT)iri^^ 
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static interactions in charged systems have been studied 
by means of mean-field approximations, such as Poisson- 
Boltzmann (PB) (or its linearized version, Debye-Hiickel) 
theory 54,55 . Despite the fact that in many instances, 
the PB theory provides a successful approximation, it has 
been shown rigorously |56U57U58| that such mean-field the- 
ories can only predict repulsion between like-charged ob- 
jects. Recent investigations have, however, elucidated the 
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important role of electrostatic correlations, neglected in 
the mean-field approximation, toward understanding the 
mechanism of like-charge attraction. Several proposals have 
been put forward in order to incorporate ionic correlations. 
These attempts include integral-equation approach J2H1 
|2"T| , perturbative improvement of the PB theory 28 29,30 
I^32I33I:MI35I36I3713SI33] . local density-functional the- 
ory 40 41 , structural-correlation approach 42 43 44 ll45| . 
harmonic-plasmon approximation |46|. and the strong-coupl 
theory 07 , 48 (that will be used in this paper) , which gen- 
erally compare well with simulation results and all exhibit 
attraction. Nonetheless, they are different in technical as- 
pects, and have their own range of rigor and applicability. 

In general, the importance of electrostatic correlations 
can be quantified by means of the coupling parameter 
S = 2Trq 3 £ 2 3 a s 07], which depends on the charge valency 
of counterions q, surface charge density of macroions <r s , 
and the Bjerrum length £b — e 2 / (AiresokBT) (associated 
with a medium of dielectric constant e and at temper- 
ature T). The PB theory is asymptotically obtained in 
the limit of vanishingly small coupling parameter S — > 
(37) . while, non-mean-field features emerge in the oppo- 
site limit of large coupling parameter, H > 1, and are ac- 
companied by strong accumulation of counterions in the 
proximity of charged objects, where counterions may also 
form strongly-correlated structures, such as a quasi-two- 
dimensional liquid or possibly a Wigner crystal at low tem- 

For rod-like and spherical macroions, effective interac- 
tions are also influenced by the entropy-driven counterion- 
condensation process, which basically arises as a result of 
the specific interplay between energetic and entropic con- 
tributions in these geometries. Clearly, a dominant attrac- 
tion between like-charged macroions necessitates a suffi- 
ciently large number of counterions being condensed in 
their close vicinity. The counterion-condensation process 
at rod-like macroions is controlled by the Manning param- 
eter [21] defined here as £ = qlsT, where r stands for the 
single-rod linear charge density. For small Manning pa- 
rameter, counterions de-condense leading to a bare elec- 
trostatic repulsion between charged rods regardless of the 
coupling parameter. While for sufficiently large values of 
Manning parameter, an attractive force may prevail as a 
certain fraction of counterions remains bound to the rods 
even in the absence of confining boundaries. An interesting 
problem is, therefore, to examine the regime of Manning 
parameters, where an effective attraction emerges. This 
has been addressed in a number of previous works for a 
pair of like-charged rods. It was argued by Ray and Man- 
ning that attraction sets in, when the Manning pa- 
rameter becomes larger than the threshold of countcrion 
condensation in the two-rod system, i.e. £ = 1/2. While, 
a different attraction range of £ > 2 was proposed by the 
counterion-condensation theory of Arenzon et al. 44 . Nu- 
merical simulations, on the other hand, reveal attraction 
between two like-charged rods in a wide range of Man- 
ning parameters including £ sw 1 (for moderate coupling 
parameters) (13U23U24U25) . though have not yet specified 
the threshold value. 



For a system of charged spheres, attraction is induced 
only when the system is confined in a box, because the en- 
ergetic binding of counterions to charged spheres is weaker 
than the entropic gain resulting in complete de-condensation 
of counterions in the absence of confining boundaries. The 
influence of the confinement volume on the effective inter- 
action between like-charged spheres was addressed in the 
simulations by Gr0nbech- Jensen et al. |14| . It was shown 
iit^at the pair potential of mean force develops a local min- 
imum at small surface-to-surface separation, which corre- 
sponds to a closely-packed bound state of spheres main- 
tained by a strong short-ranged attraction, as indicated 
independently by other simulations (15116) . This attrac- 
tion regime is separated from a repulsion regime at large 
separations by a pronounced potential barrier (14) . Upon 
increasing the confinement volume, the strength of attrac- 
tion and the barrier height decrease leading to pure repul- 
sion between spheres in a very large confining box. Other 
simulations on many-sphere systems (18U19U20"] have also 
displayed formation of large aggregates of highly-charged 
spheres, which indicate a phase separation in the ther- 
modynamic limit (1411181119) . As discussed by Linse et al. 
(19) . the equilibrium surface-to-surface separation of the 
attracting spheres in the compact state appears to be of 
the order of the counterion diameter. 

In this paper, we investigate the effective interaction 
between a pair of like-charged rods and a pair of like- 
charged spheres using the strong-coupling (SC) theory 
introduced in previous works (47) . In contrast to previ- 
ous applications of the strong-coupling theory 47 , 48 , we 
deal with curved surfaces and thus, an additional param- 
eter is introduced, which measures the ratio between the 
radius of curvature of macroions to the associated Gouy- 
Chapman length. For charged rods, this parameter is equiv- 
alent with the so-called Manning parameter (introduced 
above), and thus, we shall refer to this additional param- 
eter as Manning parameter for both charged rods and 
charged spheres. We shall identify the regime of Man- 
ning parameters, where attraction is predicted between 
macroions with or without any confinement, and exam- 
ine the equilibrium separation of interacting macroions 
for varying Manning parameter and confinement volume. 

The SC theory is obtained by a systematic virial ex- 
pansion in the limit of large coupling parameter, S — > oo, 
which is complementary to the mean-field regime (Sec- 
tion I2J • This, in particular, covers the low-temperature 
regime (but still above the crystallization temperature 
(48) ) . Nonetheless, the theory involves finite entropic con- 
tributions, thus finite-temperature effects from counteri- 
ons that appear at leading order in a simple one-particle 
form. This allows for reproducing the de-condensation pro- 
cess of counterions at low Manning parameters and there- 
fore, leads to a consistent picture in the whole range of 
Manning parameters. For sufficiently large Manning pa- 
rameters, when counterions effectively condense around 
macroions, we find a long-ranged attraction, which orig- 
inates mainly from the energetic attraction mediated by 
counterions intervening in a narrow region between macroions. 
The effective force is dominated by such an attraction 
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for increasing Manning parameter (or decreasing tempera- 
ture) and tends to a temperature-independent value, which 
scales with the macroion separation, D, as ~ 1/D for two 
charged rods, and as ~ 1/D 2 for two charged spheres. (For 
two charged plates, the effective force is obtained to be 
independent of D |37l3g] .1 In this regime, two macroions 
are predicted to form a closely-packed bound state with 
a surface-to-surface distance of the order of the coun- 
terion diameter in qualitative agreement with numerical 
findings |14U15U18II19II23II24) . Such an agreement is estab- 
lished, when the excluded- volume repulsion between coun- 
terions is negligible compared to the electrostatic interac- 
tions. For point-like counterions, the equilibrium surface- 
to-surface distance is predicted to be of the order of the 
Gouy-Chapman length -see Eqs. Ij39|l and l|58(l below- and 
a quantitative agreement is obtained between SC predic- 
tions and simulation results as reported elsewhere 25 . 

We shall also study the influence of the confinement 
due to an outer box on the effective interaction and on 
the onset of attraction. It will be shown that by decreas- 
ing the Manning parameter or increasing the size of the 
confinement, macroions undergo an unbinding transition. 
For two charged rods, the minimal Manning parameter 
to obtain attraction (the threshold Manning parameter) 
tends to the universal value of £ c = 2/3 for increasing box 
size, which is somewhat larger than the corresponding on- 
set of counterion condensation £ = 1/2. The unbinding 
transition in this system occurs continuously and is char- 
acterized by a power-law behavior. In contrast, charged 
spheres display a discontinuous transition and a potential 
barrier in the effective free energy in qualitative agreement 
with simulation results |14| . It turns out that the attrac- 
tion threshold for two like-charged spheres monotonically 
increases with the size of the confining box yielding a pure 
repulsion between unconfined spheres, but its dependence 
on the box size is logarithmically weak. Such a weak de- 
pendence of the onset of attraction on the box size may ex- 
plain the stability of aggregates of highly-charged spheres 
in large confinements that has been observed in recent 
simulations 18,19,20 . The range of applicability of the 
present theory to systems with finite coupling parameter 
will be examined by means of scaling arguments. 



2 Strong-Coupling Theory: General formalism 

Consider a system of macroions with charge distribution 
tr(r) (in units of the elementary charge e), and a num- 
ber of N oppositely-charged counterions of diameter a c 
and charge valency q. All charges interact via Coulombic 
interaction v(r) = l/|r|, and the electroneutrality condi- 
tion is assumed to hold globally in the system. Hence, one 
always has 



Nq = J dra(r). 



(1) 



Note that q and cr(r) are defined to be positive. We shall 
suppose that macroions are fixed. (In the forthcoming cal- 
culations in Sections and 01 we shall consider excluded- 
volume interactions of hard-core nature between particles- 



see Section l2~Tl l Assuming a surface charge density of a s 
for macroions, the Gouy-Chapman length, /i, is defined as 



1 



2nq£ B a-s 



(2) 



which sets a characteristic length scale for the considered 
system. Here, Ib = e 2 /(AneeoksT) is the Bjerrum length, 
the distance at which two elementary charges interact with 
thermal energy fcgT in a medium of dielectric constant e. 

It is convenient to construct a dimensionless formal- 
ism, which may be achieved by rescaling spatial coordi- 
nates, say r, with the Gouy-Chapman length as f = r//x. 
Other parameters are rescaled accordingly, for instance, 



(3) 



In the rescaled units and for point-like counterions (er c = 
0), the grand-canonical partition function of the system, 
Za, is mapped to a one-parametric field theory (apart 
from rescaled geometrical parameters such as the typical 
macroion size), and may be expressed as a functional in- 
tegral over a fluctuating field 4>(r) |4"7) 



Z A 



^exp{-7^]/~}, 



where 



S = 2irq 



3»2 



(4) 



(5) 



is the coupling parameter, Z v = yDetv contains self- 
energy contributions and the action reads 



n\6] = 



i 

2^ 



dr 



1 



(V</>(f)) 2 - uf>(i)a{T) - AQ{r)e 



(?) 



J 6 ) 

Here, A is the rescaled fugacity, and the function Q{r) 
takes geometrical constraints into consideration, e.g. re- 
stricts the positions of mobile counterions to an appropri- 
ate region in space. The electroneutrality condition, Eq. 
([T)l. written in rescaled units, relates the coupling param- 
eter to the average number of counterions through 



Q 



dra(r) = 2ir~N, 



(7) 



where Q is simply the rescaled area of macroions covered 
by electric charges. 

In general, statistical quantities such as counterions 
density profile, osmotic pressures and effective forces can 
be obtained from Eq. J3J. But, as exact calculations based 
on the partition function Q are difficult, one may consider 
limiting cases, where approximate methods are applica- 
ble. For instance, at small coupling parameters S — > 0, 
a saddle-point approximation may be applied, since the 
functional integral in this limit is dominated by the sad- 
dle point of the action, Eq. JSJ. This approximation yields 
the so-called mean-field Poisson-Boltzmann (PB) theory, 
which may be then extended to finite couplings by means 
of a systematic loop expansion around the saddle-point 
solution |37]. However, the loop expansion around PB the- 
ory was shown to be a weakly convergent series and can 
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not improve the PB results to be applicable for systems 
of large coupling parameter 0Hj • 

At large coupling parameters S 3> 1, a series expan- 
sion can be obtained for the partition function (0J) in pow- 
ers of the rescaled fugacity, A/S, as |4T) 



2 ^ = 2 »Ei 



3=0 ' 



A 
2^E 



ii i r 

n 

fc=i 



x exp{— S v(i 



dr k Q(Y k ) 



(8) 



which is nothing but a virial expansion with respect to the 
counterionic degrees of freedom. The zeroth-order term of 
the expansion, Zq, represents the partition function of the 
system, when all counterions are taken out and only the 
fixed charge distribution is present, 



-Uo/vE 



(9) 



(|4")l and JHJ, however, demonstrates that loop expansion 
and virial expansion are indeed the two asymptotic limits 
of the same problem: While the PB theory and loop ex- 
pansion become asymptotically exact in the limit of weak 
coupling S — > 0, the SC theory and the virial expansion 
produce asymptotically exact results in the complemen- 
tary limit of large coupling S — > oo. One may notice that 
higher-order terms in Eq. © contain integrations over the 
so-called Mayer function exp(— Ev) — 1, that are weighted 
by other factors coming from the interaction of counte- 
rions with macroions. Such integrations are known to be 
divergent for bulk systems such as electrolytes, where the 
virial expansion fails to produce a convergent low-density 
expansion for the equation of state. In contrast, here coun- 
terions are considered at charged macroscopic objects. It 
has been shown explicitly for charged planar walls that 
in the limit of S — > oo, all orders of the virial expan- 
sion produce finite contributions to the density profile of 
counterions and the effective force between walls IT7I. 



where 



U = ±[ dvdv'a{v)v{v-v')d{v')-^ 



2.1 The strong-coupling free energy 

dr v(r — ro)a(f ) To calculate the strong-coupling free energy, we start from 
(10) 



is the rescaled zero-particle interaction energy. (Note that 
fo is an arbitrary point chosen to fix the reference con- 
figuration with respect to which the interaction energy is 
calculated.) 

The first-order term in Eq. |(SJ| is the partition func- 
tion of the system in the presence of a single countcrion 
interacting with the fixed charge distribution via the one- 
particle interaction 



the grand-canonical free energy (in units of fcsT), 

Q A = -\nZ A , (12) 

where Z A is given by Eq. JSJ). The Legendre transforma- 
tion 

F N = NlnA + Q A , (13) 

provides us with the canonical free energy !Fn. The rescaled 
fugacity, A, in Eq. (|T3|l is calculated from 



u(r) = — — [ dr'[v(r — ?') — u(?o — r')]cr(r') 
2n J 



(11) 



N = A 



d\nZ A 
dA : 



(14) 



Higher-order terms in Eq. JHJl involve the two-particle in- 
teraction v(r — f), that emerges in the non-perturbative 
form of exp(— Sv). 

The canonical strong-coupling (SC) theory is obtained 
as an asymptotic theory from the above expansion, Eq. 
JHJ, in the limit of S — > oo. As we shall see below (Sec- 
tion 12. l| . the leading-order term in the canonical free en- 
ergy contains counterionic contributions only up to the 
one-particle terms. This reflects the fact that for large E, 
counterion-macroion correlations dominate over counterion- 
counterion correlations. Higher-order corrections for finite 
couplings have been studied analytically for systems com- 
posed of planar charged walls in Refs. [J2J021 and will not 
be considered in the present study. 

Clearly, the virial expansion, Eq. JHJ, could have been 
obtained directly from the original partition function of 
the system. Also note that although the field-theoretic 
representation of the grand-canonical partition function, 
Eq. was obtained for point-like counterions, the virial 
expansion, Eq. (|SJ), is quite general and can systematically 
include finite counterion size (a c ^ 0) via excluded- volume 
interactions between particles (see Section r2~T|) . Taking the 
detour over the field-theoretic formulation based on Eqs. 



in terms of S and other rescaled geometrical factors. (Note 
that N can be eliminated using Eq. J7J.) In general, we 
may propose the following expression for A in the large 
coupling limit S 3> 1, 



A\ A 2 
A = A + + -| 



(15) 



where Aq, A\, . . . are determined from Eq. (|14fl using Eqs. 
0, ©-(ED ^6l]. It is easy to verify that, for instance, 



A = 
Ax = 



Q 



J drfl^e-^) 

Q 2 J dfdf / ^(f)^(f / )e-"W- fi ( f ')[l 



(16) 



-Hv{r—r' 



2tt 



J drn^e-"^ 



7) 



Inserting this into Eq. l|13f) . we find the canonical free en- 
ergy Tjq (in units oiksT and up to an irrelevant additive 
constant), which also admits a large-coupling expression 



T. 



N 



^1 + ^1 



(18) 
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where the coefficient of the leading-order term is 



T\ = 7T ln 

7T Z7T 



dr J?(f )e 



-fi(r) 



+ Co, (19) 



in which Cq — (Q/2iv) InQ — Q/2ir is a constant, and the 
geometry function Q, as introduced before, specifies the 
accessible volume for counterions. 

The free energy coefficient T\ 1 which we may refer to 
as the rescaled SC free energy of the system, yields the 
actual strong- coupling free 



-r-SC _ 



(20) 




which generates the leading (non-vanishing) contribution 
to the effective forces between macroions This will 
be used in the following Sections to investigate the effec- 
tive interaction between two charged rods and two charged 
spheres in the strong-coupling limit. The SC free energy 
has a non-trivial temperature dependence that will be dis- 
cussed briefly in the Appendix for some asymptotic cases. 

For realistic systems with finite coupling parameter, 
higher-order corrections to the SC free energy become 
increasingly important as in this situation, counterion- 
counterion correlations become comparable to the counterion- 
macroion correlations accounted for by the asymptotic 
term. Thus the range of applicability of SC predictions at 
moderate to large couplings should be examined by cal- 
culating higher-order corrections. These calculations are, 
however, difficult to perform analytically for the system 
of two charged rods and two charged spheres considered 
in this paper. Thus, we shall use qualitative schemes to 
identify the regime of validity of SC results based on the 
fact that electrostatic correlations induced by counterions 
between apposing macroions are dominant, when the typ- 
ical lateral distance between counterions becomes larger 
than or comparable to the macroions surface-to-surface 
separation. This criterion, which was first introduced by 
Rouzina and Bloomfield nas been derived by calcu- 
lating higher-order terms in Eq. JSJl for planar charged 
walls 07], and is supported by numerical simulations |48) . 
This has also been addressed in simulations on charged 
rods and charged spheres in Refs. [T5lll8lll9U25| . 

Finally, as discussed before, specific excluded-volume 
interactions between particles can be accounted for sys- 
tematically within the virial expansion, Eq. (jSJ) , which al- 
lows for treating counterions of finite diameter, o~ c , within 
the strong-coupling scheme. However, as clearly seen from 
Eq. I|19|l . only macroion-counterion excluded-volume in- 
teraction enters in the leading-order contribution to the 
free energy in the SC limit. The excluded-volume inter- 
action between counterions themselves appears only in 
the higher-order terms, which are not considered here. 
For the sake of simplicity in the forthcoming calculations, 
we shall assume a hard-core repulsion between counteri- 
ons and macroions, therefore, overlapping configurations 
do not contribute to the spatial integral in Eq. <|19|l . In 
other words, the effect of counterion size in the SC limit 
is incorporated into the geometry function f], which, for 



Fig. 1. The geometry of the system of two similar rods con- 
fined in a rectangular box (top view) as considered in the text 
(Section|2J. Dashed circle shows the closest approach distance 
between macroions of radius Ro and counterions of diameter 
(T c (see Eq. i^Tf ). 



cylindrical and spherical macroions of radius i?o, implies 
a hard-core radius of 



R = R + a c /2. 



(21) 



Now, for a given amount of macroion charge, Equation 
(I21|l implies a reduced surface charge density, and thus an 
increased Gouy-Chapman length, Eq. as compared to 
the case with point-like counterions. The following results 
are, however, presented in units of the Gouy-Chapman 
length and the counterion diameter explicitly appears only 
when the actual units are restored (see e.g. Eqs. (|39|l and 
JSHJ below). 



3 Two Like-Charged Rods 

Consider two similar and infinitely-long rods each of ra- 
dius i?o, length H and the linear charge density r (in unit 
of the elementary charge e). We assume that the axes of 
the rods are parallel and located at separation D from 
each other. (We choose the frame of coordinates such that 
its origin lies in the mid-way between the rods and its 
z-axis is parallel to the rods axes - see Figure ^) To con- 
sider a more general case, we suppose that in lateral direc- 
tions, the system is confined in a rectangular box with the 
edge size of L. The electroneutrality condition is fulfilled 
globally due to the presence of oppositely-charged coun- 
terions of valency q inside the box. Assuming that the 
electric charge is distributed uniformly over the surface of 
the rods, their surface charge density reads a s — t/(2t:R), 
and the Gouy-Chapman length, Eq. @, is 



R 



(22) 



where R is the hard-core radius defined according to Eq. 
(|21Jl . As it follows from Eq. H22J1 . the rescaled rod radius, 
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R = R/ fx, is identical with the Manning parameter, i.e. 

R = f = i B qr. (23) 

The rescaled SC free energy of the two-rod system, T\, 
follows from Eq. I|19l) . in which the zero-particle interac- 
tion energy, Uq, and the one-particle interaction energy, 
u, are obtained from Eqs. (|10fl and Ijllfl respectively. Ne- 
glecting irrelevant self-energy terms and taking the off-set 
point Tq = 0, we have 



where 



^ = GnR 2 In D, 
H 

u(x,y) = 2R[]nf 1 +lnf 2 -ln(D 2 /4)}, 

h = [{x + b^f + f] 1 ' 2 , 

rz = [(x-b/2T + tf]V\ 



(24) 
(25) 



are radial distances from the rods axes. Inserting the above 
expressions into Eq. H19fl . we end up with the SC free en- 
ergy (per unit length of the rods, H, and up to an irrele- 
vant additive term) as 



where 



-J- = -2R 2 \nD- 2RlnI, 
H 



I(D,R,L) = / dxdy f2 e -2«Pafi+iiifti] 



(27) 



(28) 



The geometry function fl — fl(x, y; D, R, L) specifies the 
region of integration, i.e. it is one inside the volume sur- 
rounded by the rectangular box and the cylindrical rods, 
and is zero elsewhere. 



3.1 Threshold of attraction 

The first term in the rescaled SC free energy, Eq. 1)2 7[). 
simply gives the long-ranged (bare) electrostatic repulsion 
between the rods, and the second term involving the spa- 
tial integral contains energetic and entropic contributions 
of counterions on the leading order, which can generate 
an effective rod-rod attraction. In order to examine the 
counterionic contribution, and the onset of attraction, let 
us consider first the limit of two unconfined rods, i.e. when 
the box size becomes infinitely larger than the rod radius, 
L/R — > co. In this limit, the counterionic integral /, Eq. 
H28|) . diverges for Manning parameters R < 1/2, which can 
be seen by rescaling the spatial coordinates with the box 
size according to x — > x/L, etc., that yields / ~ L 2 ~ 4R . 
This indicates that R — 1/2 is the threshold of counterion 
condensation in this system, below which the distribution 
function of counterions around the rods vanishes |63j in 
agreement with the value obtained from Manning con- 
densation theory 3jJ . Moreover, the counterion-mediated 
force, 

F ci ~4,ln/, (29) 
oD 



vanishes for R < 1/2 as the box size tends to infinity. 
Thus, attraction between unconfined rods can only set in 
for Manning parameters larger than 1/2. Now assuming 
that the SC free energy has only one local minimum, which 
is indeed the case as we shall see from the numerical solu- 
tion, the onset of attraction can be determined by exam- 
ining the large-separation behavior of the free energy, i. e. 
for D ^> D m i n , where D m i n — 2R is the smallest possible 
axial separation of the rods. To this end, we rescale the 
spatial coordinates with the axial separation D as 



D 



y 



y_ 
b' 



Accordingly, the integral / in Eq. 128|) scales as 



I(D,R) = b 2 - 4 *J(^,R), 



(30) 



(31) 



(26) 

where J is a dimensionless integral given by 



J = 2 



x'>0 



dx'dy'[(x'+^Mm- A [(i'-^ 2 Mm-^', 

(32) 

where we have made use of the symmetry property of the 
integrand upon the reflection with respect to the plane 
x' = 0, and thus the corresponding integral runs over the 
half-space x' > excluding a disk of radius R/D centered 
at (i' = +1/2, y' = 0) (this is formally accounted for by 
the geometry function SI'). For very large D/D min , the 
radius of the disk tends to zero and the limiting behavior 
of the integral in Eq. (|32[) is determined by the contribu- 
tions from the boundary regions, which vary depending 
on whether R is smaller or larger than 1. 

For Manning parameter R < 1, the contribution of the 
boundary region around the disk vanishes, and the integral 
in Eq. I|32|l is dominated by its outer boundary, which gives 
only a constant independent of D. Therefore, the prefactor 
of J in Eq. I|31(l yields the leading D-dependence of / for 
large axial separations, and substituting this into Eq. (|27|l 
gives the limiting form of the rescaled SC free energy as 



H 



-2R(2-3R) InD. 



(33) 



For R > 1, on the other hand, the integral in Eq. (|32|) is 
dominated by the boundary region around the disk yield- 
ing J ~ £)2i?-2 f or ver y i ar g e D/D m i n , which leads to the 
following attractive tail for the rescaled SC free energy, 



^-^2R 2 \nD. 
H 



(34) 



Therefore, as clearly seen from Eqs. (|33|) and l|34|) . two 
unconfined rods experience a repulsive force at large sep- 
arations, when Manning parameter, R, becomes smaller 
than the threshold 

Rc = \, (35) 

and on the contrary, they attract each other for larger 
Manning parameters R > R c . This result is obtained also 
by numerical calculation of the free energy ( Section |3.2|l . 
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Fig. 2. The rescaled strong-coupling free energy of the two-rod 
system, Eq. I27L as a function of the rescaled axial separation, 
D — D/n, for Manning parameter R = 1.0 and different box 
sizes indicated on the graph. 



The result that the onset of attraction for two uncon- 
fincd rods is predicted to be larger than the threshold of 
the counterion condensation (R = 1/2) could be antici- 
pated from the beginning, since right at the threshold of 
the condensation, there remains an unbalanced bare re- 
pulsion between the rods according to the SC free energy, 
Eq. I|27(l . The onset of attraction has also been considered 
in a number of recent studies. Analysis of Ray and Man- 
ning 0U] based on the standard counterion-condensation 
model predicts attraction for R > 1/2. In contrast, re- 
sults from the counterion-condensation theory of Arenzon 
et al. 01] suggest a larger value of R = 2 for the threshold 
of attraction. Numerical simulations, on the other hand, 
indicate the appearance of attraction for a wide range of 
Manning parameters including R 1 13 23 : 24 , 25 . The 
precise location of the threshold of attraction between two 
charged rods is currently being studied by means of MD 
simulations (241125] . 

The attractive force predicted in the strong-coupling 
limit is long-ranged and scales inversely with the axial sep- 
aration between the rods (see Eqs. 1)33(1 and 1(34(1 and Fig- 
ureEJ- It should be noted that the SC attraction originates 
mainly from the energetic interaction mediated by coun- 
terions sandwiched between the rods, and in this respect, 
is different from the attraction obtained by the Gaussian- 
fluctuation theories 28 31 33,34 . Specifically, for increas- 
ing Manning parameter (or decreasing temperature), the 
strength of SC attraction increases and saturates to a 
temperature-independent value (see Appendix I A. 1(1 . This 
can be seen also from the large-separation behavior of the 
free energy for highly-charged rods, Eq. ((34(1 . which gives 
the attractive force as 



Frods 

H 



2iree D 



(36) 



in actual units, where F ro d s is obtained from the actual SC 
free energy, Eq. JU, as F rods = -(k B T)dF% c ' /dD. Such 



fx 
H 



-123 



-124 



i///=10 

i///=30 

i/ju=100 

» v L/fi= oo 



20 



40 



D/fx 



60 



80 



100 



Fig. 3. The same as Figure but for Manning parameter 
£ = 0.3. 



an energetic attraction is also obtained for planar charged 
walls |47lf4"5| and charged spheres (Section^) , but with dif- 
ferent dependencies on D. The mechanism of SC attrac- 
tion for large R, therefore, qualitatively agrees with the 
low-temperature results |44U45| (see also Refs. (34U38U46) 
for discussions on the crossover from the low-temperature 
to the high-temperature regime). 



3.2 Equilibrium axial distance 

Now, let us consider a more general case with a rectan- 
gular box of edge size L (Figure 0. For sufficiently large 
Manning parameters R > R c = 2/3, the presence of con- 
fining walls is expected to have only a minute effect on 
the effective rod-rod interaction since most of the counte- 
rions are strongly localized in the proximity of the rods. 
Figure [3 shows the rescaled SC free energy of the two-rod 
system, calculated numerically from Eq. I(27ll . as a func- 
tion of the axial separation, D, for Manning parameter 
R = 1 and various box sizes. As seen, the free energy 
quickly converges to the free energy of an unconfincd sys- 
tem with L = oo, and the long-ranged attraction is only 
slightly strengthened at large separations, when rods are 
located close to the box boundaries. The free energy takes 
its minimum at a small axial distance, -D* ~ 2i?, which 
corresponds to a closely-packed bound state of rods. (We 
may refer to as rescaled "equilibrium" axial separa- 
tion.) The approximate form of T\ in the vicinity of its 
minimum can be calculated as 



3\ 
H 



6R 2 \n D - 2R\n(D ~ 2R) 



(37) 



for R 3> 1 and up to an irrelevant additive term (Appendix 
IA.1() . The first term in Eq. ((37(1 contributes a dominant 
(energetic) attractive force and the second term generates 
a repulsive force between the rods at small separations. 
Using Eq. ((37|) . the rescaled equilibrium axial separation 
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Fig. 4. Rescaled equilibrium axial separation of two like- 
charged rods plotted as a function of Manning parameter (the 
rescaled rod radius), R = R/fi, for the box sizes L — L/fj, — 
10, 100, 1000, oo as indicated on the graph. The small arrow 
locates the threshold Manning parameter 7? c =2/3, and the 
thick grey curve corresponds to the contact axial separation 



is obtained approximately as 



•J K 



(38) 



for sufficiently large Manning parameter. Equation l|38|l 
shows that the (rescaled) equilibrium surface-to-surface 
distance of the rods, £>* — 2R, decreases and tends to a 
value of the order of unity as Manning parameter tends to 
infinity. Restoring the actual units in Eq. (|38[1 . we obtain 
£>* « 2R + 2^/3 + 0(n 2 /R), where R and /j, are given by 
Eqs. (|21(1 and (|22(1 . Thus, the actual equilibrium surface- 
to-surface separation is obtained approximately as 



73* — 2Rq 



(39) 



when the Manning parameter is sufficiently large (or the 
Gouy-Chapman length, /i, is small), that is about the 
counterion diameter. 

For small Manning parameters R < R c = 2/3, qual- 
itatively different features arise. The effective interaction 
between the rods exhibits a significant dependence upon 
the size of the confining box and the location of the min- 
imum of the free energy tends to infinity for L — > oo re- 
sulting in a pure repulsion between unconfined rods (see 
Figure . This occurs as a result of the dilution of the 
counterion cloud around weakly-charged rods. For R< 1, 
the equilibrium axial separation scales with the box size 
and is obtained approximately as (Appendix IA.1|) 



L 

\/7T 



(40) 



Now, for a given box size, the free energy varies smoothly 
from the typical forms shown in Figure to those shown 
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Fig. 5. Rescaled equilibrium axial separation of two uncon- 
fined rods (L — oo) in the vicinity of the threshold Manning 
parameter R c = 2/3 (solid curve with symbols). The dashed 
line shows a power-law fit with the slope a = 3/2 (see Eq. 

ED)- 



in Figure [3] by decreasing Manning parameter. There is, 
however, a rapid change near R c = 2/3 indicative of a 
continuous unbinding transition in an infinitely large box. 
This is seen more clearly in Figure 01 where the predic- 
tions of the SC theory for the equilibrium axial separation 
is shown as a function of Manning parameter for various 
box sizes. (These results have been obtained directly by 
numerical minimization of the free energy, Eq. Ij27(l .) The 
region above a given curve corresponds to the axial separa- 
tions at which the rods attract, and the region below that 
corresponds to the repulsion regime. When the box size 
tends to infinity, the equilibrium axial separation diverges 
for R — > i?+ revealing the power-law behavior 

D*~(R-R c )- a (41) 
shown in Figure where within our numerical errors 

a = 3/2. (42) 



3.3 The regime of applicability of SC results 

The system of two like-charged rods has been investigated 
in recent numerical simulations revealing strong attrac- 
tive forces at moderate to large Manning parameters |13l 
I23H24U25| . The SC theory provides quantitative predic- 
tions for the equilibrium axial separation, that may be 
compared with simulations results. However, such com- 
parison with realistic systems of finite coupling parame- 
ter (e.g. S - 10 - 100 for simulations in Refs. [T3II23II24I 
125) ) must be performed in a certain regime of parame- 
ters, where the present asymptotic theory (strictly valid 
only for 5 — * oo) is applicable. To this aim, extensive 
simulations have been performed recently by Arnold and 
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Holm |24| covering a wide range of parameters, which al- 
low for examining the influence of higher-order electro- 
static and excluded-volume effects. The comparison with 
these results will be presented in separate publications |24l 
I25| , and in the following, we merely introduce two criteria 
identifying the regime of parameters, where a reasonable 
agreement between simulations and SC predictions is ob- 
tained. 

For highly-charged rods, excluded-volume interactions 
between counterions may be significant as attraction is 
accompanied by accumulation of counterions in a narrow 
region between the rods fAppcndix lA.lfl . In this situa- 
tion, the typical distance between intervening counterions 
(lining up in the z-direction) , a z , may be estimated from 
the local electroneutrality condition, q = 2ra z , giving 



q_ 

2r' 



(43) 



where q is the counterion valency and r is the single-rod 
linear charge density. When, the counterion diameter, cr c , 
is larger than a z , 

a z <(T c , (44) 

the excluded-volume repulsions between counterions be- 
come important, and further accumulation of counterions 
between the rods according to the SC mechanism is pro- 
hibited. As a result, the equilibrium separation of the rods 
(determined from the electrostatic component of the total 
force between them) appears to be larger than the present 
predictions, since the strength of the mediated attraction 
drops. Though, it should be noted that volume interac- 
tions may indeed contribute an additional attractive com- 
ponent to the total force for highly-charged systems [221 
I24| . It is a priori not clear which of these effects is stronger. 

On the other hand, as frequently quoted in the lit- 
erature |16lFT 8 19 42 . 47 . 48 . electrostatic correlations in- 
duced by counterions between two apposing macroions are 
dominant, when the typical distance between counterions, 
a z , becomes larger than the macroions surface-to-surface 
distance, A = D — 2Rq, i.e. when 



A < a z 



(45) 



where a z is given by Eq. In rescaled units, this con- 
dition reads 



where the coupling parameter, Eq. J2J), is 



-,3s2 



R 



(46) 



(47) 



4 Two Like-Charged Spheres 

Consider a system of two like-charged spheres each of ra- 
dius i?o and charge valency Z. The spheres are located 
at a center-to-center separation of D, and together with 
their neutralizing counterions are confined in a cubic box 
of edge size L. (The geometry of the system is similar to 
what we have sketched in Figure^ where the two solid cir- 
cles now display the largest cross-section of the spheres.) 
We assume the electric charge to be distributed uniformly 
over the surface of the spheres, thus their surface charge 
density reads u s = Z/(4irR 2 ), and the Gouy-Chapman 
length, Eq. J3J), is 

9 R 2 

(48) 

where R is the (hard-core) radius of spheres, Eq. J2J. In 
rescaled units, we have 



R 



R 
/' 



£ B qZ 
2R ' 



(49) 



which is referred to in the following as Manning parameter 
for charged spheres in analogy with the cylindrical case, 
Eq. J22I). 

To study the strong-coupling interaction of the two 
spheres, we follow similar lines as presented in the preced- 
ing Section by calculating the rescaled SC free energy, Eq. 
(|19f) . The zero-particle and one-particle interaction ener- 
gies are obtained from Eqs. illOfl and (|llf> respectively as 



R 

U {) = -28tt-^ 
D 

u(i,y,S) = ~2R 2 [-^ + ^-Z- 



r 2 D 



(50) 



where 



h = [{Z + D/2) 2 + ?+~z 2 ] 1 l 2 , 
r 2 = [{x-D/2f+f+~z 2 ]V 2 , 



(51) 



are radial distances from the centers of the spheres. The 
reference point is arbitrarily chosen to be in the mid-way 
between the spheres i"o = 0, and self-energy terms are 
neglected. Using Eqs. i[T9"j) and (JSDJl, the rescaled SC free 
energy is obtained as 



d4 

Tx = -AR 2 lnI, 
D 



(52) 



where 



for charged rods. The above criterion identifies the regime 
of parameters, where basically the counterion-rod corre- 
lations are superior to the counterion-counterion correla- 
tions and the present asymptotic theory is applicable |251 
I47U48| . For moderate coupling parameters S ~ 10, this 
regime can be reached using charged rods with Manning 
parameter larger than £ ss 1 |25j . 



I(D,R,L) 



dxdydz f2 e 



(53) 



The geometry function fi(x, y, z; D, R, L) specifies the re- 
gion accessible for counterions, i.e. it is one inside the 
volume surrounded by the cubic box and the two spheres, 
and is zero elsewhere. 
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The first term in Eq. I|52l) gives the long-ranged bare 
repulsion between the spheres, and the second term in- 
volving the spatial integral, /, contains energetic and en- 
tropic contributions of counterions on the leading order 
reproducing the counterion-condensation process and an 
attractive force for highly-charged spheres. 

It follows easily from Eqs. (|52|l and (|53[l that for fi- 
nite Manning parameter, the counterion-mediated force 
vanishes as the size of the confining box tends to infinity 
resulting in pure repulsion between unconfincd spheres. 
This goes back to the fact that in the absence of confining 
boundaries, spherical macroions can not bind counterions 
that tend to diffuse in order to gain entropy. This is re- 
flected by the divergency of the counterionic integral, Eq. 
<|53|) . as L — > oo. (Note that the involved integrand is al- 
ways positive and larger than one, and upon rescaling the 
spatial coordinates as f = r/L, it scales with the volume 
of the confining box as / ~ Z 3 .) As a result, the distribu- 
tion function of counterions 63 as well as the counterion- 
mediated force, 

F ci ~4,m/, (54) 
oD 

that scales with the box size as ~ L~ 2 , vanish for increas- 
ing box size. In contrast, it turns out that two like-charged 
spheres can attract each other in any finite box provided 



that their Manning parameter is sufficiently large as will 
be demonstrated below. 

Figures E^L-d show the typical form of the SC free en- 
ergy, calculated numerically (using Monte-Carlo integra- 
tion methods) from Eq. (|52|l , for two spheres in a box with 
L = 100 and for several Manning parameters. As seen, 
at sufficiently small Manning parameter, the long-ranged 
repulsion is dominant and only a weak attraction oper- 
ates, when the spheres are at large separations comparable 
with the box size (FigureEK)- But at Manning parameters 
larger than a threshold Rc (which is obtained as Rc ~ 
3.3 for L — 100), a local minimum (short-ranged attrac- 
tion) emerges at small separations about D w 2R corre- 
sponding to a meta-stable bound state of spheres (Figure 
EJd). The attractive interaction develops further and the 
potential barrier disappears, when the Manning parame- 

~ (2) 

ter becomes larger than the second threshold R c (which 

is R { c ] « 3.8 for L = 100) - see Figure Hi. This clearly 
demonstrates the occurrence of a discontinuous transition 
between a closely-packed bound state and a repulsion- 
dominated state of two spheres by changing Manning pa- 
rameter R. 

In Figured we show the locations of the minima (solid 
curves) and the maximum (dashed curve) of the SC free 
energy as a function of Manning parameter for L = 100 
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Fig. 8. The threshold Manning parameter plotted as 

a function of the box size L (see Eq. I|590 in the text). This 
threshold is obtained for box sizes larger than L « 20. 



Fig. 7. Locations of the minima (solid curves) and the maxi- 
mum (dashed curve) of the SC free energy for the two-sphere 
system plotted as a function of Manning parameter R for 
L = 100. The two arrows on the horizontal axis show the loca- 
tions of the threshold Manning parameters R c and R,i ■ The 
dot-dashed line represents the contact separation D m in = 2i?, 
and the vertical dotted line shows the Manning parameter for 
which the values of the two minima of the free energy are equal. 



(obtained numerically using Eq. Q52JI). The upper branch 
represents the shallow repulsion-dominated minimum, which 
is strongly sensitive to the box size. For small Manning 
parameter i!« 1, the corresponding equilibrium center- 
to-center distance increases linearly with the box size and 
is obtained approximately as 



which yields an equilibrium surface-to-surface distance, 
— 2R, of the order of unity in rescaled units. Now 
restoring the actual units in Eq. (|57J) . we obtain £>* w 
2R + An/7 + 0(^ 2 /R), where R and H are given by Eqs. 
(I21|l and (|48|l . Thus, the actual equilibrium surface-to- 
surface separation of spheres is obtained approximately 
as 

A* = D* - 2R Q w <j c + -h + 0(li 2 ), (58) 



7 



for L 



when the Manning parameter is sufficiently large (or the 
Gouy-Chapman length, /i, is small), that is about the 
counterion diameter. 

Similar features are obtained for box sizes that are 
larger than roughly L ss 20; for smaller box size only the 
shallow large-distance minimum is obtained as depicted in 
Figure 01. The striking result is that the threshold Man- 
oo (AppendixE21). By contrast, the small-separation^ng parameters R { c 1} and R {2) monotonically increase by 



(55) 



minimum (the lower branch) is effectively independent of 
the box size and is maintained essentially by an attrac- 
tive interaction mediated by counterions intervening be- 
tween the spheres. For sufficiently large Manning param- 
eter R ^> 1, the approximate form of the rescaled SC 
free energy around its minimum at small separations is 
obtained as 



increasing the confinement volume. This behavior is shown 
in Figure [8] for R^ (symbols), along with the best loga- 
rithmic fit to these results (solid curve) as 



RP =a + blnL, 



(59) 



- 28 * 

D 



AR 2 hi(D - 2R) 



(56) 



up to an unimportant additive term (Appendix IA.2|) . For 
increasing Manning parameter (or decreasing tempera- 
ture), the free energy expression fBfUl leads to an (ener- 
getic) attractive force independent of temperature, which 
scales with the center-to-center distance of the spheres, 
D, as ~ 1/D 2 , in agreement with results obtained by 
Shklovskii 05]. Using Eq. (|56|l . the rescaled equilibrium 
center-to-center distance is obtained approximately as 



D*^2R 



; + o(±>, 



(57) 



where a « 0.55 and b w 0.6. Such a weak increase of the at- 
traction threshold demonstrates that even in a very large 
confinement, two like-charged spheres can fall into the at- 
traction regime provided that the corresponding Manning 
parameter exceeds a moderate threshold. This behavior 
may also explain the stability of closely-packed clusters 
of highly-charged spheres in large confinement, and their 
insensitivity to the box size as addressed in recent sim- 
ulations |14II16II18II19| . The simulation parameters, where 
strong attraction has been observed, indeed covers the pre- 
dicted attraction-dominated regime R > R { c 1] (see Table 

One should note that the unbinding transition can oc- 
cur also by changing the confinement volume for a given 
Manning parameter. This is implied by the fact that, for 
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Simulation 
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Co (A) 


L(A) 
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R 




S(i) 


ax (A) 


A s im (A) 


A, (A) 


Gr0nbech- Jensen 


2 


10 


7.01 


7 


3.3 


50-200 


1.07 


8.1 


26 


2.8-3.7 


7.7 


2.5 


3.41 


et al. [U 
Wu et al. (T5| 


2 


20 


7.14 


10 


4 


100 


1.01 


11.9 


28 


3.3 


7.5 


4 


4.10 


Allahyarov 


2 


32 


112 


48.9 


4.4 


~ 10 2 


0.73 


70.1 


615 


~ 3.3 


25.5 




4.41 


et al. [JJ] 

Linse et al. |18I19| 


3 


60 


7.15 


20 


4 


~ 10 2 


0.75 


29.2 


85 


~ 3.3 


9.8 


4 


4.05 


Hribar et al. jUj 


3 


12 


7.15 


10 


2 


~ 10 2 


0.94 


11.7 


68 


- 3.3 


11.0 




2.09 



Table 1. Parameters from simulations on highly-charged spheres (see the text for the definitions). The last two columns show 
the equilibrium surface-to-surface distance, A s i m , obtained in these simulations (if explicitly measured), and the strong-coupling 
prediction, A, — D, — 2Rq. Some of the numbers are given up to the order of magnitude, and the extracted values of A from 
simulations have a typical resolution of about lA. 
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Fig. 9. Rescaled equilibrium center-to-center separation of Fig. 10. Rescaled equilibrium surface-to-surface separation of 

like-charged spheres as a function of Manning parameter in the like-charged spheres as a function of Manning parameter in the 

attraction-dominated regime (i.e. R > Ri^). Symbols show attraction-dominated regime. Symbols show data from recent 

data obtained from recent simulations (circle: Gr0nbech- Jensen simulations (circle: Gr0nbech - Jensen et al. [JJ|, square: Wu et 

et al. [TI], square: Wu et al. \Ej\, triangle: Linse et al. |18I19| 1 al. PJ, triangle: Linse et al. |18ll9p and the solid line is the 

compared with the SC prediction (solid line). SC prediction. Dot-dashed line shows the contact separation 



a given Manning parameter between the two thresholds 

FtP < R < R [ c \ the height of the free energy barrier 
decreases slowly by increasing the box size. Qualitatively 
similar features have been obtained in the simulations by 
Gr0nbech- Jensen et al. ^1] on a system of two like-charged 
spheres with divalent counterions. The obtained potential 
of mean force between spheres was shown to have a local 
minimum at small separation (the compact state) sepa- 
rated by a potential barrier from a repulsion regime (see 
Fig. 1 in Ref. 14 , which schematically resembles Figure|SJ 
Note that the tail of the effective interaction is influenced 
by the boundary conditions, that is periodic in the case of 
the simulations). The barrier height was found to decrease 
for increasing confinement volume resulting in repulsion 
between spheres (Fig. 3 in Ref. 14 ). The appearance of a 
potential barrier between two highly-coupled spheres was 
also noted in Ref. |21|. 



The strong-coupling regime has also been investigated 
in several other simulations both on many-sphere systems 
pilT9ll2T1j and on a single pair of spheres [TSITH] . Tabled 
presents typical parameters from some of recent simula- 
tions. In the two last columns of the Table, we compare 
the equilibrium surface-to-surface separation between at- 
tracting spheres from these simulations, A s i m (if explic- 
itly estimated), with the SC prediction, A* = — 2i? , 
in actual units. (In Figures El and we also show the 
center-to-center distance, and the surface-to-surface 
distance, D r — 2R= (£>» — 2R — <r c ) / in rescaled units.) 
The strong-coupling results have been calculated numeri- 
cally from the free energy, Eq. I|52|l . and involve the finite 
size of counterions, which as discussed in Section l2.ll en- 
ters at the leading order via a finite closest approach dis- 
tance between spheres and counterions. As seen, there is a 
qualitative agreement between the theoretical predictions 
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and the simulations results. In fact, as addressed in pre- 
vious numerical studies (see for example the Discussion 
in Ref. |19j). the equilibrium surface-to-surface separation 
of spheres in these simulations appears to be of the order 
of the counterion diameter. This result agrees with the 
SC prediction for highly-charged spheres of small Gouy- 
Chapman length (jj, ~ lA) as demonstrated by the limit- 
ing expression 1)58)1. 

The regime of applicability of SC predictions to sys- 
tems with finite coupling strength, such as the above sim- 
ulations, may be examined by considering the ratio of 
the lateral counterion separation at charged spheres, a±, 
to the surface-to-surface distance between the spheres, 
A = D - 2R |16II18II19II42II47II48| . The typical value of 
a± may be estimated from the local electroneutrality con- 
dition ira 2 ^ ~ q/cr s , which holds up to a geometric factor 
of the order of one. In rescaled units, we get aj_ / /i = dj_ ~ 
(23) 1 ' 2 , where the coupling parameter, Eq. (JjjJ, reads 



for charged spheres. The lateral separation between coun- 
terions may also be written as 




As discussed before, the SC regime, i.e. when higher-order 
corrections are small (Section HJ, is characterized by the 
condition 

A < a±. (62) 

The estimated values of aj_ for the given simulations are 
shown in Tabled As seen, the above criterion is fulfilled 
for the simulation parameters [1611181119) . In addition, the 
size of counterions are typically smaller than their lat- 
eral separation indicating that excluded-volume repulsion 
between counterions, which enters through higher-order 
corrections to the SC theory, is not a dominant effect. 

5 Conclusion and Discussion 

In this paper, we have presented the predictions of the 
strong-coupling theory for the effective interaction between 
two like-charged rods and two like-charged spheres. The 
SC theory obtained from a systematic virial expansion 
(Section^ provides a framework to study strongly-coupled 
systems (3 3> 1). For large Manning parameters (or weakly- 
curved macroion surfaces), counterions effectively condense 
around macroions and a long-ranged attraction is obtained. 
(For charged rods, attraction is obtained even in an in- 
finitely large confining box due to a universal counterion- 
condensation effect in this limit, while spheres show at- 
traction only when they are confined in an outer box.) 
In general, the SC attraction is accompanied by strong 
accumulation of counterions in the intervening region be- 
tween macroions that gives rise to attractive forces, which 
scale with the macroions separation, D, as ~ 1/D for two 
rods, and as ~ 1/D 2 for two spheres. The equilibrium 



surface-to-surface separation of highly-charged macroions 
is then predicted to be of the order of the counterion di- 
ameter plus a term of the order of the Gouy-Chapman 
length -see Eqs. ()39)1 and 1)58)1 - which is in qualitative 
agreement with results obtained from recent simulations 
on both charged rods [21] and charged spheres (Section^). 
For point-like counterions, a quantitative agreement has 
been established with recent Molecular Dynamics simu- 
lations [5SJ. Some of the features such as the finite-size 
dependence of the effective interaction and the appear- 
ance of a potential barrier in a system of charged spheres 
are also consistent with the simulation results [14U21) . 

Note that in the simulations, attractive interactions 
are found to have a short-ranged character, whereas the 
SC theory predicts long-ranged attractions for large Man- 
ning parameters. This difference should be attributed to 
finite-coupling effects as in the simulations considered here 

only moderately large, 3 ~ 10 — 10 2 . As it was demon- 
strated for planar charged walls |47U48| . for such coupling 
parameters, the strength of electrostatic correlations is 
substantially decreased when charged walls are far from 
each other and consequently, the range of attraction is 
decreased. (In that case, it was shown that the large- 
distance behavior of the effective interaction is described 
by the PB theory and is purely repulsive). Yet, the SC 
results were shown to remain valid at sufficiently small 
separations, as obtained from the criterion explained in 
the text [23114 71148) (see Eqs. gSJ and ftEfy). As we al- 
ready discussed, this criterion can be satisfied at moder- 
ate to large couplings using moderate values of Manning 
parameter, e.g. R > 1 for two charged rods |25) . Exper- 
imentally, this regime can be reached using multivalent 
counterions and highly-charged macroions. For instance, 
in aqueous solutions of DNA (with radius i?o ~ 10 A and 
linear charge density re ~ 6e/nm), Manning parameter 
and the coupling parameter are as large as R ~ 8 and 
3 ss 25 in the presence of divalent counterions, and as 
R « 12 and 3 m 80 for trivalent counterions such as sper- 
midine. In colloidal dispersions, an aqueous solution of 
highly-charged surfactant micelles of, for example, typical 
radius Ro w 20A and charge valency Z 60 represents 
Manning parameter and the coupling parameter of the or- 
der of R « 30 and 3 « 100 for trivalent counterions. It 
is also noted that recent experimental observation of at- 
traction between like-charged colloidal particles in Refs. 
64 : 65 66 concerns the regime of small coupling parame- 
ters (typically with 3 ~ 10~ 2 — 10 _1 due to the large size 
of spheres ~ 1/im). These observations, therefore, should 
not be compared with the strong-coupling results. 

We have also examined the predictions of the SC the- 
ory for small Manning parameter (or small radius of curva- 
ture), e.g. for the onset of attraction. For sufficiently small 
Manning parameters and increasing confinement volume, 
one expects that counterions completely de-condense lead- 
ing to pure repulsion between macroions. It should be 
noted, however, that when counterions de-condense, counterion- 
macroion correlations become effectively small and contri- 
butions from higher-order corrections to the asymptotic 
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theory may become important even for large coupling pa- 
rameters (this can be understood also in terms of the crite- 
rion introduced in Sections l3.3l andl4^l. These effects might 
result in finite corrections for the predicted threshold of 
attraction or scaling exponents, such as Eq. JI2J|, in the 
limit of large couplings. However, as we showed, the de- 
condensation process at small Manning parameters is cap- 
tured by the SC free energy and a dominant repulsive in- 
teraction is obtained, which gives a consistent picture for 
the whole range of Manning parameters. Note also that 
the threshold of condensation obtained for two charged 
rods (Section EH} coincides with the value obtained from 
Manning theory |49ll59j . Further numerical studies in this 
regime will be useful to investigate the role of higher-order 
corrections. (These aspects have been studied for a single- 
rod system by means of Monte-Carlo simulations |60|.l 

Finally, we emphasize that the thermodynamic behav- 
ior of macroionic solutions and colloidal dispersions may 
not be derived directly from the SC free energy presented 
in Sections[3]and01 As we showed, by decreasing Manning 
parameter or increasing the box size, attracting macroions 
undergo an unbinding transition to a repulsion-dominated 
state, which occurs continuously for two rods and discon- 
tinuously for two spheres. These transitions do not nec- 
essarily represent thermodynamic phase transitions, be- 
cause in our study, the distance between macroions is as- 
sumed to be fixed, while in realistic situations (as well as 
in the simulations on charged spheres cited in Section ^J, 
the macroion-macroion separation is an annealed degree of 
freedom contributing separately to the partition function. 
Also, for fluctuating spheres there will be a logarithmic 
contribution to the entropy for large separations due to 
the increasing free volume available for the sphere-sphere 
distance coordinate. On the other hand, there has been ev- 
idence of a thermodynamic phase transition of the first or- 
der (phase separation) from recent numerical simulations 
on highly-charged spheres [TIHT^HT Q , 20 . The systematic 
study of such phase transitions in the strong-coupling limit 
will be subject of a future investigation. 
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A Asymptotic analysis of the SC free energy 

A.l The two-rod system 

Let us first consider the limit of vanishingly small Manning 
parameter R <C 1. In this case, the integral in Eq. I|27(l may 
be expanded in terms of R, and the rescaled SC free energy 
per unit length, H, of the rods is subsequently obtained 
as 

5- « - 2R 2 ln£) + ^2- ( 1 dxdy [lnfj + lnf 2 ] 

H L 2 J-L/2 



- 2R\nL 2 , (63) 

where we have omitted terms that are independent of 
L and D as they are irrelevant in our discussion. The 
asymptotic expression (|63|l involves the rescaled mean en- 
ergy of the system of two-rods and neutralizing counteri- 
ons (the first two terms) together with the entropic con- 
tribution of the counterfoils (the last term), which has 
an ideal-gas form. (This is seen more clearly using Eqs. 
(I2U[1 and (|63|l and restoring the actual units that yields 
Tfp « -21 b t 2 H \nD + 2£_N(\n n + In r 2 ) - N In L 2 per 
ksT and up to some additive terms. Note that Man- 
ning parameter equals the rescaled radius of rods, Eq. 
(I23|l . thus R — > implies line charges in the rescaled 
picture.) This form of the SC free energy shows that at 
small Manning parameters, the counterions are effectively 
unbound, though still interacting with the charged rods. 
The bare repulsive force between the rods coming from the 
first term scales like ~ D^ 1 , but the countcrion-induced 
attraction force coming from the second term scales as 
~ D/L 2 , which becomes vanishingly small as L — > oo. 
Therefore, the asymptotic free energy l|63l) admits a shal- 
low L-dependent minimum as it is seen in Figure By 
minimizing expression l|63|) with respect to D, the bound- 
state separation is approximately found as in Eq. I|4U|) for 
L — > oo. 

For large R > R c , as seen in Figure |3 the free energy 
substantially decreases, when the rods are close to each 
other. Inspection shows that in this situation, the main 
contribution to the integral I in Eq. (|27|l comes from the 
intervening region between the rods (-D/2 + R < x < 
D/2 — R; y w 0). In fact, this is associated with strong 
accumulation of counterions in this region for increasing 
Manning parameter, as it can be checked directly from 
the counterionic distribution function obtained in the SC 
limit The integral /, Eq. (J2SJ, may be rewritten as 

I = J dxdyf2e- 2ka ^-y\ (64) 

where 

g(x,y) = [In f i + In r 2 ], (65) 

and f\ t 2 are defined in Eq. I|26|l . It turns out that (x,y) = 
(0, 0) is the saddle point of g, thus for sufficiently large R, 
we may use a saddle-point approximation to calculate J, 
which gives (up to some irrelevant prefactors) 

jRJe -a«mzj> x £2 x / dxe 8Ri \ (66) 

J-1/2+R/D 

where we have rescaled the coordinates as x — > x/D, y — ► 
y/D (Eq. i|30|0 and assumed that the box size is suffi- 
ciently large. Indeed, the above approximation is valid 
only for sufficiently large R/D ratio, since g is singular 
at (x/D, y/D) = (±1/2,0). In other words, it remains 
valid as long as the surface-to-surface separation of rods 
is sufficiently small, i.e. 5 — (D — 2R)/2R <C 1. In this sit- 
uation, the following approximate expression is obtained 
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for the rescaled free energy using Eqs. Ij66() and (|27Jl . 

Si . . n 62 1 



re 6R z lnD -2Rf(D,R), 



(67) 



where / ps ln(D — 2i?) + £)(<5). Now, by minimizing expres- 
sion i|67|) with respect to D, the equilibrium axial separa- 
tion is obtained as in Eq. I|38l) for R 3> 1. We emphasize 
that Eq. (|67|l represents the approximate form of the free 
energy around its local minimum, and is not valid for large 
separations (or small Manning parameters close to the at- 
traction threshold R c — 2/3). The large-distance behavior 
of the free energy is given by Eqs. I|33|l and l|34|l . 

Restoring the actual units in Eq. (|67|l . and noting that 
the rescaled free energy T\ is related to the actual SC free 
energy (per k B T), Fff, through Eq. we find 



Fn G ~ 6£ B r 2 HlnD - N\n(D - 2R). 



(68) 



The first term in this equation formally corresponds to 
the energetic attraction between two rods mediated by 
neutralizing counterions that are located between them 
(lining up in z direction with an effective linear charge 
density of 2r). While the second term may be regarded 
as an (repulsive) entropic contribution from counterions. 
Note that the attractive component of the resultant force 
is independent of the temperature, which only appears 
as a prefactor in the second term. Clearly, for increas- 
ing Manning parameter (or decreasing temperature), the 
strength of attraction increases |13II34[ and saturates to 
the maximal value given by the first term in Eq. (|68(l . 



A. 2 The two-sphere system 

For vanishingly small Manning parameter R <C 1, one 
may expand the integral involved in the rescaled SC free 
energy, Eq. (|52|l . which yields 



-4 



8R 4 f L/2 - A 
/ d 6 r\- 
L 3 J-L/2 ri 



i] 

r 2 



AR 2 In L 3 . 



(69) 



This expression is basically similar to the low-Manning- 
parameter expansion for the two-rod system, Eq. (jHSJ, 
containing the mean energy of the two spheres with neu- 
tralizing counterions (the first two terms) together with 
the entropic contribution of the counterions (the last term) 
that has an ideal-gas form. (This is seen using Eqs. i|20[l 
and Q69[) and restoring the actual units that gives Tfj C ~ 
£ B Z 2 /D - N£ B Zq(l/n + l/r 2 )-N In L 3 in units olk B T.) 
The bare repulsive force coming from the first term scales 
like ~ D~ 2 , but the counterion-induced attraction force 
coming from the second term scales as ~ D/L 3 , which be- 
comes vanishingly small as L — > oo. The asymptotic free 
energy (|69|l has a shallow L-dependent minimum as seen 
in Figures Eli-c, the location of which is approximately 
obtained as in Eq. (|55|l in the text. 



For large Manning parameter, as seen in Figure|S[i, the 
free energy exhibits a local minimum at small center-to- 
center separations, where the main contribution to the vol- 
ume integral in Eq. 152|) comes from the intervening region 
between spheres (-D/2 + R < x < j)/2-B;y« 0; z ps 0). 
In this saddle-point approximation similar to the 

two- rod system (Eq. I|66|) ) may be performed, which gives 
the following asymptotic expression for the rescaled free 
energy for sufficiently large R and small surface-to-surface 
separation 6 = (D — 2R)/2R < 1, 



R 

Fx = -28^ -AR 2 f[D,R), 



(70) 



where / ~ \n(D — 2R) + 0(S). The equilibrium separation 
is subsequently obtained by minimizing expression (|70|l 
with respect to D yielding Eq. I|57|l in the text. 

Using the free energy expression (|70|l and restoring the 
actual units, the approximate form of the actual SC free 
energy, Eq. I)2()[l. is obtained as 



sc 



Z 2 



\ - -7£ B — -Nln(D-2R). 



(71) 



The first term in Eq. I|71|l formally corresponds to the 
energetic attraction between two spheres mediated by a 
neutralizing counterion (of charge valency q = 2Z) located 
between them. The second terms may be regarded as an 
(repulsive) entropic contribution from counterions. Upon 
increasing Manning parameter (or decreasing the temper- 
ature), the strength of attraction increases and saturates 
to the maximal value given by the first term in Eq. I|71() . 
which gives a temperature-independent force that scales 
with the center-to-center distance as ~ 1/D 2 . 
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